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ABSTRACT 



Aims. We determine the number density and area contribution of small-scale inter-granular Call bright G-band structures in images 
of the quiet Sun as tracers of kilo-Gauss magnetic flux-concentrations. 

Methods. In a 149" x 117" G-band image of the disk center at the activity minimum, 7593 small inter-granular structures were seg- 
mented with the 'multiple-level tracking' pattern recognition algorithm. The scatterplot of the continuum versus the G-band brightness 
shows the known magnetic and non-magnetic branches. These branches are largely disentangled by applying an intrinsic CallH ex- 
cess criterion. The thus obtained 2995 structures contain 1152 G-band bright points (BP) and 1843 G-band faint points (FP). They 
show a tendency toward increasing size with decreasing G-band excess, as expected from the 'hot wall' picture. Their CallH and 
G-band brightness are slightly related, resembling the known relation of Ca II and magnetic field strength. The magnetic flux density 
of each individual BP and FP is estimated from their G-band brightness according to MHD model calculations. 
Results. The entity of BP and FP covers the total FOV with a number density of 0.32/ Mm^ and a total area contribution of 2.0%. 
Their individual calibrations yield a mean flux density of 20 Mx/cm- in the entire FOV and 13 Mx/cm^ for inter-network regions. 

Key words. Sun: magnetic fields - Sun: activity - Sun: photosphere. 



1. Introduction 

The small-scale magnetic flux-concentrations on the non-active 
solar surface are an important ingredient for understanding so- 
lar magnetism. Active region features, such as sunspots, faculae, 
prominences, and flares, exhibit a strong variation with the solar 
cycle, commonly explained by a global solar dynamo acting in 
the deep convection zone. On the other hand, small-scale mag- 
netic flux-concentrations exist that cover the solar surface more 
homogeneously in time and space indicating a less tight depen- 
dence on the global dynamo 

Unfortunately, it is difficult to measure magnetic flux on 
those scales predicted by MHD modeling (Vogler & Schiissler 
2007). Measurements of circular polarization due to the Zeeman 
effect can only give the net flux, since opposite polarities cancel 
each other out within the spatial resolution element (Dommguez 
Cerdena et al. 2006). The absolute (unsigned) flux, however, has 
to be deduced from magnetically broadened intensity profiles 
(Stokes-I) that are afifected not only by the magnetic field but 
also by the, a priori unknown, intrinsic atmosphere of each indi- 
vidual flux concentration. 

The problems for such a conversion have been discussed by, 
e.g., Sanchez Almeida (2003). Algorithms used for this proce- 
dure are often called 'inversion codes', which might suggest 
that the observed polarization profiles are actually 'inverted' into 
magnetic field parameters. Instead, these codes vary both the 
magnetic field structure and the intrinsic atmosphere model until 
an optimal parameter fit between calculated and observed pro- 
files is found. The codes also optimize the filling factor of the 
magnetic signal inside the resolution element but do not con- 
sider two-dimensional radiative effects in very small structures. 
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Recent observations impressively show that flux- 
concentrations can be of such a small size that even the 
0.09" spatial resolution of the 1 m SST at La Palma (Scharmer 
et al. 2003) does not reach the lower end of the size distribution 
of G-band structures. Stokes images, however, cannot achieve 
such a high spatial resolution, due to additional optical elements 
and longer exposure times. Hence, they even depend more on 
the (doubtful) fill factor than intensity images. 

Many of the problems in interpreting polarization measure- 
ments may be avoided by analyzing intensity images (Sanchez 
Almeida et al. 2004; De Wijn et al. 2005). Small-scale magnetic 
concentrations of kilo-Gauss flux density are known to exhibit 
excess brightness in the CH-band near 430 nm (G-band) over the 
neighboring continuum (Berger et al. 1995), thus allowing their 
study in post factum reconstructed images of higher spatial and 
temporal resolution than achieved for Stokes measurements. The 
excess in the G-band brightness over the neighboring spectral 
continuum is a different measure than the excess over the mean 
photospheric G-band level, defining G-band bright points (BP; 
Muller & Roudier 1984). The latter cover two families of struc- 
tures: (I) inter-granular BP, which are magnetic (Berger & Title 
2001) and down-drafting (Langhans et al. 2002), and (II) non- 
magnetic BP, which move upwards with the respective granules 
(same two papers). 

We separate both types of BP using a pattern recognition 
code, (Bovelet & Wiehr 2007) especially adapted to detect the 
inter-granular BP. We then select magnetic inter-granular struc- 
tures using their G-band-to-continuum ratio and additionally 
their CallH excess brightness. There is some controversy about 
these criteria: Rezaei et al. (2007) find that the Call excess can- 
not be considered a general indicator for magnetic flux concen- 
trations. Berger et al. (1998) argue that the G-band brightness is 
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Fig. 1. Magnetic map of the solar disk on April 12, 2007 from 
MDI showing an extremely low level of solar activity. The 
149" X 1 17" FOV investigated in this paper is marked. 



not a sufficient criterion for the magnetic nature of a BP. De Wijn 
et al. (2005) suggest the Call brightness to be a more realis- 
tic criterion. Their concern about the G-band criterion refers to 
Berger & Title (2001), who, however, only exclude that criterion 
for those BP that appear embedded in granules, i.e., the above 
family II. 

We avoid these, restricting our selection to small struc- 
tures, thus excluding extended granular features and, hence, non- 
magnetic local G-band enhancements embedded. Among the 
small-scale inter-granular structures (IgS), we remove those that 
populate the non magnetic branch in the G-band-to-continuum 
scatterplot. The remaining IgS are than sampled by a significant 
Ca II excess. 

The powerful pattern recognition algorithm, followed by 
a selective sampling of small structures, elimination of non- 
magnetic ones by their G-band-to-continuum ratio, and classi- 
fication of magnetic ones by the Ca II excess, qualifies our study 
as an alternate way to investigate small-scale solar kilo-Gauss 
magnetic flux-concentrations at the highest possible spatial res- 
olution and independent of interpreting polarimetric measure- 
ments. 



2. Segmentation of G-band features 

We analyzed a set of high-resolution images, observed at an ex- 
ceptionally deep solar activity minimum on April 12, 2007 (cf., 
MDI-magnetogram Fig.[TJ. The Dutch Open Solar Telescope 
(DOT) was pointed at the intersection of the central meridian 
with the solar equator (0" E/W; 6" N) by P. Sutterlin, who took 
image bursts simultaneously in the G-band, CallH, Ho-, blue 
and red continuum, which he subjected to image reconstruction 
following de Boer (1996) and finally composed to a 4-image mo- 
saic of 149" X 1 17" total field-of-view. That large FOV, marked 
in Figure[T]and shown in Figure|2]as an Ha image, covers 1/221 
of the total solar surface, allowing a statistically significant in- 
vestigation of small-scale magnetic flux-concentrations. 




Fig. 2. Reconstructed Ha image of a 149" x 117" region cen- 
tered on the central meridian and the equator (6''N) on April 12, 
2007, (blue and red wing images added; courtesy P. Sutterlin, 
Utrecht). The square marks a subregion used for a visualization 
of the data processing in Fig.|5] 



As a first step, the G-band image was analyzed using the 
MLT_4 recognition codeQ (Bovelet & Wiehr 2007) to obtain an 
appropriate pattern of solar structures and their individual sizes 
in pixel counts. The latter depend to some extent on the 'uni- 
tary cut level' equally applied to each of the segmented cells: 
the higher the cut level, the smaller the size. Our final choice 
of a unitary cut level of 0.5 was adapted to a realistic represen- 
tation of IgS sizes, accepting slightly broadened inter-granular 
lanes. As a second step, we selected inter-granular structures 
('IgS') with an area < 75 pixels (2 x lO'^cm^). This procedure 
excludes larger features with local G-band enhancements occa- 
sionally embedded that are known to be non magnetic (Berger & 
Title 2001; Langhans et al. 2002). The suitability of the achieved 
representation is demonstrated in Figure|5] for a sub-field indi- 
cated in Figure|2] 



3. Selection of IgS with Ca ii H excess 

For the thus selected total of 7593 IgS, we plot the G-band ex- 
cess over the intensity in the neighboring spectral continuum fol- 
lowing Berger et al. (1998). The spatial correlation between the 
two simultaneously observed G-band and (blue) continuum im- 
ages is optimized by a cross-correlation to one pixel (0.071"). 
Figure[3] shows a scatterplot of the 33% highest intensity pixels 
of each IgS in the continuum and in the G-band. The 33% restric- 
tion (rather than its mean value) ensures that the deduced IgS 
brightness is largely independent of individual shapes. Figure[3] 
gives both brightness values in units of the mean photosphere at 
FOV locations free of local enhancements (e.g., the lower-left 
region in Fig.|2]). 

The scatterplot in Figure[3] shows two branches that, ac- 
cording to measurements by Berger et al. (1998) and to calcu- 
lations by Shelyag et al. (2004), represent magnetic and non- 
magnetic IgS and overlap for G-band brightness /'-''' < 1 . 1 x Q„„ 
(shaded in Fig. [3]). Since we do not, a priori, know which of the 
segmented features are magnetic, we remove, as a third step, 
3150 non-magnetic IgS in the upper part of the non magnetic 
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branch (sample A), located above the kinked line line inserted in 
FigureO 
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Fig. 3. Scatterplot of blue-continuum and G-band brightness for 
all 7593 IgS selected, grouping in a magnetic (lower) and a 
non-magnetic branch that are separated by the full line; overlap 
regime shaded. 



Among the remaining 4443 IgS (sample B), we disentangle 
non-magnetic and magnetic ones by their Ca II H excess bright- 
ness, assuming magnetic IgS to exhibit a substantial CallH 
excess, whereas non-magnetic IgS are considered to have no 
significant CallH brightening (Lites et al. 1999; Rezaei et al. 
2007). We correspondingly remove, as a fourth step, the non 
magnetic IgS from sample B applying a threshold of the intrinsic 
Ca n H excess defined as 



ExCa 



with I up denoting the CallH brightness contrast over the mean 
for 33% of the brightest and //o„. for 33% of the darkest pixels of 
each IgS. Application of the ratio of both contrasts as the expo- 
nent to lup avoids excluding features with a high but flat intrin- 
sic Ca II H brightness distribution, as G-band bright points. Such 
a criterion mainly selects IgS with a spatially isolated CallH 
excess inside their corresponding G-band pattern and excludes 
faint and 'fuzzy' CallH structures (widely considered to have 
an acoustic nature). 

The influence of this intrinsic Ca II H excess criterion on the 
number density of IgS is demonstrated in Figure|4] For sam- 
ple A, the logarithmic number density of IgS already drops at 
very low ExCa values, since this group is mainly populated with 
IgS from the non magnetic branch without intrinsic Ca II H ex- 
cess. They mark granular fragments or faint and 'fuzzy' CallH 
brightenings that are already removed by a low ExCa value, 
since both basis and exponent of ExCa are small. 

For sample B, the logarithmic number density of IgS remains 
constant up to ExCa x 0.9, then drops, and finally linearly de- 
creases for ExCa > 1 .4. The high Ca II H excess of IgS in sam- 
ple B causes already low ExCa thresholds to exclude most of 
the non-magnetic IgS, then leaving the number density constant. 
When further increasing the threshold to ExCa > 0.9, magnetic 
IgS with Ca II H excess are also excluded progressively. 

We eventually eliminate from sample B those IgS that do 
not reach ExCa - 1.2, and obtain 2995 remaining IgS, which 
are Ca II H bright and thus reasonably assumed to be magnetic, 
hence, denoted as 'MIgS'. Among these, G-band bright MIgS 
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Fig. 4. Number density (logarithmic) as a function of the intrin- 
sic Ca II H excess for all 7593 inter-granular structures, the 3 150 
above and the 4443 below the division line in Fig[3j flat top and 
steep decrease approximated by lines intersecting at ExCa= 1.1. 



with I'^^ > 1.1 X I^,gg„ (i.e. beyond the overlap regime) are the 
known BP. Accordingly, we denote the remaining MIgS with 



1.1 X /; 



Cb 



as 'G-band faint points' (FP). 



4. Results 

4.1. Call brightening 

Figure|5] gives the visual impression that the Ca II H brightness 
is not tightly related to the G-band brightness. This is quantita- 
tively seen in the scatterplots in Figure|6l Our 2995 MIgS show a 
fairly good relation that resembles the one found by Rezaei et al. 
(2007) between the Ca II K index and the magnetic flux density. 
The smooth transition between BP and FP (at Z"^* = 1.1 x /,^*„„) 
indicates their similar (magnetic) nature, as found by Bovelet & 
Wiehr (2007). The absence of such a relation for the 3150 IgS 
from sample A (right panel) indicates that these are mainly not 
(or weak) magnetic. 

The findings by Lites et al. (1999) and by Rezaei et al. (2007) 
that magnetic flux concentrations in the network cell interior do 
not exhibit a Ca II excess might suggest that most of our MIgS 
are located at network boundaries. To check this idea, we super- 
pose the Ha image from Figure|2]with the location of our finally 
selected MIgS (Fig.|7]). This superposition impressively shows 
that the BP are preferably located near the dark Ha fibrils, which 
are widely assumed to have their roots in small-scale kilo-Gauss 
magnetic flux-concentrations. 

In contrast, FP appear also outside the Ha fibrils, i.e., in the 
inter-network (IN). This, however, means that at least part of the 
IN flux-concentrations do exhibit marked CallH excess, since 
the FP have been equally selected by the ExCa threshold. The 
discrepancy of this finding to other observations might be due to 
the higher spatial resolution achieved in our data. This is read- 
ily seen in the lowest panel of Figure|5] which displays the same 
region as the upper panels but after artificial degradation of the 
Call image simulating a spatial resolution of 0.75" (15 cm tele- 
scope aperture). It impressively shows, that even at such mod- 
erate degradation, the isolated tiny Ca II H bright features are no 
longer detectable. Our intrinsic Ca II excess criterion even fails 
to select IgS for a faint degradation to a 0.5" spatial resolution. 
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Fig. 6. Scatterplot of the brightness in Ca II H and in the G-band for the 2995 Ca II H bright IgS from sample B (MIgS) and for the 
3 150 IgS from sample A (right panel), showing a tight relation only for the MIgS (here linearly fitted); the lower limit of the ordinate 
values marks the influence of the Ca II H threshold applied to sample B . 



Hence, very high spatial resolution is required to check whether 
small-scale magnetic flux-concentrations actually exhibit a Ca II 
excess. 



4.2. Size 



The spatial resolution achieved in our data is proven by the 
power spectrum (Fig. [8]), which drops to twice the noise level 
at 2n-k « 37 Mm"', i.e., 170km on the Sun. Accordingly, the 
histograms of BP and FP in Figure|9^ drop towards small sizes, 
reaching 170km (0.23") at 5% of their maxima. The most fre- 
quent sizes of the 11 52 BP range between 1 2 and 20 pixels, those 
of the 1843 FP between 14 and 23. Since our analysis does not 
need any assumption about feature shapes, we only give diam- 
eters for equivalent circular shapes additionally at the upper ab- 
scissa of Figure|9^ and obtain most frequent circular diameters 
of (220 + 30) km for BP and (250 + 30) km for FP 

Towards increasing sizes, the histogram of BP establishes 
the existence of an upper size limit for the BP (cf., Wiehr, 
Bovelet & Hirzberger 2004; Puschman & Wiehr 2006; Beck et 
al. 2007). The corresponding diameter of an equivalent circular 
shape amounts to 400km, in agreement with Muller (1983) and 
with Bovelet & Wiehr (2003). The FP seem to be systematically 
larger We verify this by plotting the G-band brightness excess of 
the MIgS over the neighboring spectral continuum as a function 
of the size and find that larger MIgS are primarily faint (Fig.^). 

This is a hint that for larger magnetic flux-concentrations the 
'hot wall' is wider, thus reducing the brightness excess. In turn, 
small FP do exist that are of similar size as the BP (lower left part 
of the data cloud in Figure|9]3); their 'closer hot walls' must then 
be fainter, probably due to a lower magnetic flux density (cf., 
4.4). The few data points in the upper right of Figure|9]3 represent 
14 elongated MIgS (among 2995) that have not been segregated 
enough by MLT_4; this proves the power of our algorithm. 



4.3. Area contribution 

We reasonably assume that the G-band appearance of each MIgS 
is a realistic measure of the corresponding kilo-Gauss magnetic 
feature, following Shelyag et al. (2004), who state that the "spa- 
tial G-band pattern matches the spatial structure of the mag- 
netic feature". Hence, the sum of all 2995 MIgS sizes may be 
considered as an actual area contribution in the 9.2 x 10' tm^ 
FOV (3.48 X 10^ px). The 2995 MIgS cover 1.85 x 10'^ cm^ 
(697 19 px) corresponding to 2.0% of the FOV and a number den- 
sity of 0.32 Mm"^ (Table 1), i.e., the same as found by Sanchez 
Almeida et al. (2004). 

The location of BP and FP in the Ha image (Fig.|7]) allows to 
distinguish their affiliation to fibril and non-fibril (IN) regions. 
The 'voids of fibrils' appear to be homogeneously populated by 
FP, whereas the BP primarily occur in centers of fibrils regions, 
as is expected from their nature as magnetic footpoints. For a 
quantitative analysis, we approximated those complementary re- 
gions by a binary mask, obtained from sufficient smoothing of 
the Ha image. We find a number density of 0.22 Mm"^ in the 
fibril voids, i.e., about 2/3 of the 0.32 Mm"- found for the entire 
FOV (Table 1). We verified that this value does not significantly 
depend on the smoothing for a binary mask. 

4.4. Estimate of magnetic parameters 

To obtain the mean magnetic flux density of the MIgS in the 
FOV, we might multiply their area contribution of 2% by a 'typ- 
ical' value for an intrinsic flux density of 1 kG and obtain a mean 
(unsigned) flux density of 20 G. A possible way estimating the 
intrinsic flux of each MIgS may be obtained from the relation 
between the G-band brightness and the magnetic flux density, as 
modeled by Shelyag et al. (2004). There, in Fig. 4, most of the 
scatter of the G-band brightness arises from the spatially varying 
continuum level in their calculation area (also containing gran- 
ules). Our sample of MIgS, however, is restricted to small inter- 
granular features with a certain G-band excess and will be less 
influenced by continuum variations. We may thus consider the 
mean MHD relation as calibration. 
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Fig. 10. Calibration curve for MIgS (left) averaged from Fig. 4 in Shelyag et al. (2004); scatterplot of G-band brightness and 
magnetic flux (middle) and flux histogram (right panel) of 2995 MIgS. 



The modeled relation refers to each pixel of a magnetic fea- 
ture. We averaged the G-band brightness of the pixels of each 
MIgS and used that mean for calibration (left panel of Fig.fTOt. 
thus obtaining the intrinsic magnetic flux density. Multiplying 
this value with the corresponding area (returned by MLT_4) 
yields a flux value for each MIgS. The middle panel of FigurefTOl 
shows these flux values as a function of the mean G-band bright- 
ness for each of our 2995 MIgS. 

The flux values depend to some extent on the unitary cut 
level (0.5; cf.. Section 2) applied to each segmented cell after 
normalization to its brightness maximum in phase 4 of MLT_4 
(cf., Bovelet and Wiehr 2007): the higher the cut level, the 
smaller the size. Choosing a lower cut threshold would expand 
MIgS to their peripheral pixels that are fainter and thus would 
lower the mean G-band brightness of the features. The calibra- 
tion curve (left panel of Fig.fTOt will then assign lower flux den- 
sities, and the flux of this expanded MIgS will increase by only 
a few percent. The influence of our unitary cut level on the de- 
duced MIgS flux is thus negligible. 

4.5. Mean flux in the FOV 

The histogram of the flux of the MIgS (right panel of Fig.fTOt 
shows a narrow distribution with most frequent values of 0.45 x 
10'^ Mx and a total flux of 1.86 x lO^' Mx for all 2995 MIgS. If 
we divide the latter value by the total area of the 2995 MIgS 
(1.85 X lO'^cm^), we obtain a most frequent MIgS flux den- 
sity of 1005Mx/cm^, showing that the kilo-Gauss assumption 
is a realistic mean over the entire FOV. But the most frequent 
flux density largely differs between BP and FP, amounting to 
1597Mx/cm2 for the BP and to 775 Mx/cm- for the FP 

The difference between BP and FP is not that great, how- 
ever, for the mean flux density in the entire FOV. Dividing the 
total flux of the 1152 BP and that of the 1842 FP separately 
by the entire FOV area, we obtain 9.6Mx/cm^ for the BP and 
10.6Mx/cm^ for the FP (which are more frequent but less mag- 
netic, see above). The sum of both values gives a mean flux 
density over the FOV of 20.3 Mx/cm^ for all 2995 MIgS, a 
number that ranges between that by Daras-Papamargaritis & 
Koutchmy (1983; 10^' Mx/lO-^cm^) and the one by Sanchez 
Almeida (2003; 30Mx/cm2 for his fuU FOV). 



Table 1. Characteristics of IgS with ExCa > 1.2, distinguished 
for BP (/c* > 1.1 X FP (/«' < 1.1 X and MIgS 

(BP plus FP) in the entire FOV; MIgS also in regions without 
pronounced Ha fibrils. 



sampled 


total 


no. dens. 


area cover. 


mean flux dens. 


entity 


no. 


[Mm-2] 


[%] 


[Mx/cm-] 


BP 


1152 


0.12 


0.60 


9.6 


FP 


1843 


0.20 


1.40 


10.7 


MIgS in FOV 


2995 


0.32 


2.00 


20.3 


MIgS in voids 


888 


0.22 


1.52 


12.8 



In the 'fibril voids' seen in Figure|2] we find 0.22MIgS/Mm^ 
covering 1.52% of the corresponding void area (Table 1). 
Assuming kilo-Gauss field strengths for each MIgS, we obtain 
a mean flux density in the fibril voids of 15.2Mx/cm^, in per- 
fect agreement with the finding by Sanchez Almeida (2003; for 
the network interior). However, attributing to each MIgS its in- 
dividual magnetic flux via the calibration curve (left panel of 
Fig. [Tot yields a mean flux density of only 12.8Mx/cm^ for the 
void areas. This value is lower than the above one for kilo-Gauss 
fields, since the fibril voids are primarily populated by FP, which 
exhibit lower most frequent flux densities. We verified that the 
reduction of the mean flux density in fibril voids as compared 
to that in the entire FOV does not significantly depend on the 
binary mask used to exclude fibril regions (cf.. Table 1). 



5. Conclusions 

The segmentation and selection of small-scale inter-granular 
structures with a significant Ca II H excess allows lower limits 
to be determined for the number density and fractional area cov- 
erage of inter-granular kilo-Gauss magnetic flux-concentrations. 
A magnetic calibration of each feature yields a mean (unsigned) 
magnetic flux density in a quiet solar region close to existing 
results. Our method offers an alternate approach to polarime- 
try, although it inevitably misses those magnetic structures that 
do not exhibit a Call H excess and 'elusive magnetic structures' 
from fluctuations in thermodynamic properties (which Sanchez 
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Fig. 5. Sub-region of 10.4" x 20.7" marked in Fig.|2]and rotated 
by -90" in the light of the G-band, together with the correspond- 
ing patterns of bright points (BP; white), faint points (FP; light 
gray), and granules (dark gray); the CallH image without and 
with artificial degradation to 0.75" spatial resolution; tickmarks: 
IMm. 



Almeida, 2000, estimate to carry at least half of the solar mag- 
netic flux). 

In case the total flux in our disk center FOV were represen- 
tative of the entire solar surface (221 x FOV), the Sun would be 
coveredby>0.6xlO''MIgS with a total flux of >4.1 x 10^^ Mx. 
Considering also magnetic structures which are not covered by 
our selection, the Sun at activity minimum may well be covered 
by a few miUions of MIgS with a total flux of 10^"* Mx or even 
more. 




Fig. 7. Location of magnetic inter-granular structures MIgS 
within the Ha structures of the 149"xl 17" FOV; white dots mark 
G-band bright MIgS (BP), black dots mark faint ones (FP); the 
dots represent locations rather than feature sizes; gray peripheral 
zones not covered by the Ha mosaic; tickmarks: 10 Mm. 
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Fig. 8. Power spectrum of the G-band image showing a limit- 
ing resolution of 2jt x k = 37/Mm of the DOT mosaic image 
observed and reconstructed by P. Siitterlin. 



Our method offers a quantitative investigation of highly re- 
solved (reconstructed) intensity images, avoiding the ambigu- 
ities of the Stokes calibration procedure and allowing closely 
neighboring opposite polarities to be disentangled up to the spa- 
tial resolution achieved in two-dimensional brightness images. 
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age reconstruction, Dr J. Hirzberger the power spectra. An anonymous referee 
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Fig. 9. Size histograms of 1152 BP and 1843 FP from sam- 
ple B. The lower abscissa gives reaHstic sizes in pixel counts; 
the upper abscissa gives equivalent diameters of corresponding 
circular shapes. The data points represent counts per bin-size of 
21 maximum-scaled histograms with different bin-sizes that are 
overplotted to provide enough invariance from subjective sam- 
pling. The lower panel shows the G-band brightness excess (with 
respect to the neighboring spectral continuum) as a function of 
the size for all 2995 MIgS (BP plus FP). 
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